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ABSTRACT 


Helicopters are occasionally used to carry loads hanging below them 
on cables. Piloting the helicopter under these conditions is difficult, 
particularly when the mass of the load is comparable to that of the heli- 
copter, and there are gusty winds. An autopilot logic is designed here 
for controlling a helicopter with a hanging load. A 16th order model for 
the system is decoupled into four subsystems: (a) A second order system 

for yawing motion, (b) a second order system for vertical motion, (c) a 
sixth order system for longitudinal motion, and (d) a sixth order system 
for lateral motion. A novel measuring scheme, which could be used in 
remote areas, is developed and filters are designed to estimate the state 
variables from these measurements. 

The autopilot can be used to move the load over short distances 
without retracting the cables. This is done by automatically shifting 
the autopilot modes from position-hold (hover) to acceleration-hold to 
velocity-hold (cruise) to deceleration-hold to velocity-hold (near hover) 
to position-hold (hover). Use of such an autopilot might save consider- 
able turnaround time. 

The Sikorsky S-61 helicopter is chosen as an example vehicle. The 
performance of the controlled system is studied in the presence of lon- 
gitudinal and lateral winds. Satisfactory response is obtained under 
design conditions and also with nominal changes in system parameters. 
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Chapter I 


INTRODUCTION 


A helicopter is an ideal vehicle for transporting heavy, bulky 
loads over short distances where surface transport is either infeasible 
or uneconomical, e.g. , in off-loading containers from ships or carrying 
transmission towers or prefabricated buildings to remote sites. In 
many cases the load is so big that it cannot be carried inside the air- 
craft. Therefore it must be transported hanging from cables fastened 
to the aircraft. It might be carried in a hanging position to save 
loading and unloading time. 

The pendulum modes of the hanging load couple into the motions of 
the helicopter to form an unstable system. The instability is not 
severe and an experienced pilot can still control the helicopter. 
However, the pilot work load is so great that it is almost impossible to 
add further tasks such as position-hold (precision hover), even in still 
air. On the other hand, an automatic control system can stabilize the 
system and perform additional tasks such as precision-hover even in 
gusty winds. 

To simplify the autopilot design and implementation the complete 
motion is approximated by four uncoupled motions. A technique is 
developed for finding filter and controller gains. The Sikorsky S-61 
is taken as the example helicopter. 

For this system some changes in parameters are expected during a 
mission. Also some parameters may not be known accurately in advance. 
The effect of the parameters on system behavior is studied using a 
fixed set of filter and controller gains. 

Chapter II develops the equations of motion for the system and 
describes the measuring technique near hover. The system definition 
matrices for the linearized system are given under design conditions. 

Chapter III describes the design technique for controller and 
filter for the precision hover autopilot. The performance is studied 
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under design conditions. 

Chapter IV suggests a technique for moving the hanging load over 
short distances. The load transfer is carried out in four steps involv- 
ing speeding up, holding the speed, slowing down and precision hover. 

Chapter V makes a study of the effects of changes in system param- 
eters. Chapter VI summarizes the problem and the solution. 
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Chapter II 


EQUATIONS OF MOTION 


2. 1 Introduction 

A helicopter carrying a hanging load can be modeled as a system of 
three connected bodies: (a) Rotor, (b) Fuselage, and (c) Hanging load. 

The rotor can be tilted much faster than the fuselage or the cable carry- 
ing the hanging load. Since we are not interested in high frequency 
effects associated with rotor states, it suffices to use a quasi-steady 
rotor model. In this description it is assumed that the inclination of 
the rotor no-feathering plane (NFP) to the fuselage can be changed 
"instantaneously" using the cyclic pitch control. The hanging load is 
modeled as a point mass. Thus the dynamical system is simplified to one 
rigid body and a point mass (the vehicle and the hanging load).* 

This mathematical model is of 16th order with four control variables. 
The state variables are: Three position and three translational velocity 

coordinates of the vehicle center of mass; three angular orientation and 
three angular velocity coordinates of the fuselage; two angular orienta- 
tion and two angular velocity coordinates of the cable. The control 
variables are the longitudinal and lateral inclinations of the rotor 
NFP, the collective pitch, and the tail rotor. However, near hover, the 
yaw motion and the vertical motion are very nearly uncoupled from the 
longitudinal and the lateral motions. This results in a second order 
model for yaw motion with tail rotor as control, a second order model 
for vertical motion with collective pitch as control, and a 12th order 
model for longitudinal and lateral motions with longitudinal and lateral 
cyclic pitch as controls. We will later discuss the decoupling of the 
longitudinal motions from the lateral motions. 


In case the size of the hanging load is comparable to the length 
of the cable, the hanging load and the cable could be considered as a 
compound pendulum. The analysis is very similar and the number of state 
variables is the same. 


3 




4 


FIG. 2.1 Helicopter Carrying a Hanging Load 


2.2 


Equations of Motion 


Figure 2.1 shows the coordinate system used to describe the motions 
of a helicopter carrying a hanging load. It is assumed that the load is 
carried by several cables attached to the helicopter at a point below 
its center of mass. Some important distances are shown in the figure. 

Let 

x,y,h = longitudinal, lateral, and vertical deviations from a desired 
point 

u,v,w = three components of linear velocity 
<S>,6,ty = roll, pitch, and yaw angles 
p,q,r = roll, pitch, and yaw rates 

R,v,£ = subscripts for rotor, vehicle, and hanging load (or cable) 
respectively 

V ,V ,V = three components of wind velocity 
w w w 
x y z 

d. . = aerodynamic damping coefficients 
( • ) 


The kinetic energy of the system is 


1 r . 2 . 2 • 2-, 1 r • 2 • 2 • 2-, l r *2 *2 i 2-> 

T = 7 j m [x +y +h^ + — m [x.+y „+h 1 + ■=[ I 4> +1 0 +1 W ] 

2 v v J v v J 2 r £ £ £ 2 U x v y v z v 


and the potential energy is 

V - m gh + m gh 
v & v r £ 


( 2 . 1 ) 


( 2 . 2 ) 


The vertical position of the hanging load can be expressed in 
terms of other variables, thus 


h„ = h - b cos Q cos <t> - £ cos cos it> 

£ v v v £ £ 


\ - b - 1 * 1 ^ + 1 


for small angles. 
Also, 


x 

y 


£ 

£ 


x + b sin 6 + £ sin 0, a x + b0 + £6 

v v £ v v £ 


b<b 

v 




(2.3) 
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giving 


6 = (x -x -bS )/Z 

Z Z v v 

= ( “ y i + V b4> v )/i 


( 2 . 4 ) 

( 2 . 5 ) 


Thus the potential energy is 


V = (m v +ra^)gh v - m^g(b+i) 


V 


_ _ i/x-x -b0 \ / -y .fy \ 


-y.+y -b<j> 
z v v 


2. 


( 2 . 6 ) 


The independent coordinates are h y ; ; 9^ , x v > x g and 

<b . y , y Using the generalized Lagrange method, the equations 

v * J v Z 

of motion can be written. The equation governing h y is: 

m h + m h + (m +m„)g = T - m d. (h -V ) 
vv Z v v Z v h v w 

z 

Since the nominal value of vertical acceleration is zero 


T = (m +m )g 
n v Z 


( 2 . 7 ) 


where subscript n denotes nominal value. 
If 


T = T + 6 T 
n 


( 2 . 8 ) 


and 


Y = 


m +m. 
v Z 

m 


( 2 . 9 ) 


we have 



( 2 . 10 ) 


The equation governing is: 


if = T - T - I 
z v t m z W v 


( 2 . 11 ) 


where T is the torque produced by the tail rotor at the center of 

mass and T is the torque of the main rotor 
m 
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( 2 . 12 ) 



The equations governing 0 , x , and x „ (longitudinal motions) 

v v i 


are: 


mg 

10 + m gb0 + — — (x -x -b0 )(-b) 

yv £, v i, £ v v 


= Ta0 - a.I § + CT. I (x -V ) + ol I $ + CT I (y -V ) 

R 1 y v lyvw x ' T “3 y v 3 y y v 


+ Di e R + d 2 <d r 


(2.13) 


m i g 

m x + — — (x -x -b0 )(-l) = -T(0 +0 O ) + a o m 0 - CT m (x -V ) 

vv i £ v v vR 2vv 2 vvw 

x 

- a.rn i - CT m (y -V ) - D 0 - D f (2.14) 

4vv 4 v v w 3R 4 R 

y 


m 


n^g 


x. + — — [x -x -b0 ] = -m CT _ (x -V ) 
l £ £ v v J jf 5 i w 


(2.15) 


The first term in the generalized force expressions on the right-hand 
side of equations (2.13) and (2.14) is due to rotor thrust, the next 
four terms represent damping and the last two terms arise from blade 
coning, offset distance of flapping hinge, and other effects. Expres- 
sions for , D 2 , Dg , and have been derived by Hall (HI). The 

only term on the right-hand side of equation (2.15) is due to air drag 
on the hanging load. 

For small changes we can approximate T by T in the equations 

n 

above. Thus, using (2.7) and (2.9), we have: 
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The equations for 0 y , and y, (lateral motion) can be 

v v S, 

written in a similar manner: 



2. 3 Decoupling into Subsystems 

The complete system can be divided into three subsystems. They are: 
(a) Yawing Motion: This is a second order system in which l/r and 

are state variables and is the control variable. 
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d 

dt 



1 


V 

+ 

r— 

n 0 

d . 


V 


1/1 

z 


(2.18) 


(b) Vertical Motion: This is a second order system with h^ and h^ 

as state variables and 6 t , related to the collective pitch, as 
the control variable. 



0 

0 



+ 


0 

1 

Ym 


6 T 



(2.19) 


(c) 


Longitudinal and Lateral Motions: These motions are described by a 
12th order mathematical model with two controls. The state vector, 
x , is 


x = (9 ,B , x , x ,x ,x ,$ ,<j> ,y ,y ,y.,y.) , 
— v v v v i £ v v v v £ £ 


the control vector u = (6*,<J> ) and the noise vector 

- R R 

w = (V , V ) . The state equation is 
— w w 

x y 


x = Fx + Gu + Tw 


( 2 . 20 ) 


Matrices F(12tX12) , G(12x2) , and r(12X2) are given in Table II-l. 


The models describing the vertical motion and yawing motion are 
simple second order systems. It is quite straightforward to design 
control laws for these systems. We will focus our attention on design- 
ing a control law for the system describing the lateral and longitudi- 
nal motions. 


2.4 


Decoupling the Longitudinal from the Lateral Motions 


Referring to Table II-l,.we see that the equations governing the 
longitudinal and lateral motions are coupled by cross-damping torques 
and forces. If we assume that 
motions and 


6 primarily controls the longitudinal 
R 


<J> primarily controls the lateral mations, there is an 
R 

additional coupling in the closed loop system because 4> produces a 

R 
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System Definition Matrices 
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small pitching moment and a small longitudinal force (in addition to 

producing large rolling moment and lateral force), while 9 produces 

R 

a small rolling moment and a small lateral force (in addition to pro- 
ducing large pitching moment and longitudinal force). However, a 

certain linear combination of 9 and <t> produces large pitching 

R R 

moment and longitudinal force, a small rolling moment and no lateral 
force. Another linear combination produces large rolling moment and 
lateral force, a small pitching moment and no longitudinal force. If 
we take these two linear combinations as control variables, the longi- 
tudinal-lateral decoupling approximation is quite good for the purpose 
of designing controllers and filters. 


2.5 Measurements for Estimating the State Variables 

All of the state variables in the system presented in the last 
section are observable from measurements of the position of the vehicle 
center of mass (or the sling load). However, with these measurements 
alone, the errors in the angular orientation estimates are large even 
if a relatively accurate position measurement is available. Hence, it 
is almost essential to also have measurements of vehicle and cable 
angles. 


The position measurement technique suggested here consists of a 
second cable (a "measurement cable") fastened to a point on the ground, 
which passes through a ring mounted elastically in the plane of the 
aircraft floor, and is held by a constant tension winch inside the 
aircraft (Figure 2.2a), The portion of the cable above the aircraft 
floor is parallel to a reference axis fixed in the aircraft. The 
forces in the springs are measured using strain gauges, while a poten- 
tiometer in the constant tension winch measures the length of the cable. 
The nomenclature for the measurement cable is shown in Figure 2.2b. Let 

0 and 0 be the inclinations of the rope to the aircraft reference 
m m 

axis and L m , the length of the cable. Then neglecting the small 
deflections of the springs 


x i+lf - b £ = x l + y j - h £ 

m mmvv v v 


( 2 . 21 ) 


\ 
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O 
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FIG. 2. 2a Measurement Technique 
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FIG. 2.2b Nomenclature for Measurement Cable 


where 


— > 

L = 
m 


in 


V l+tan^3) +tan 2 <l> 


[-tan© , tan$ , -l] 


m 


m 


m 


m 


and the transformation between 1 , j , ft , and i » J* v » and 6 
small angular inclinations of the aircraft is 



_ 

1 t -0 


r^1 

1 

V 


V V 



J 


-t 1 <t> 


A 

J 

V 


< 

< 



£ 


0 -<t> 1 


£ 

V 


V V 




Using this transformation in equation (2.22), 


L = 
m 


m 


Vl+tan 2 ® +tan 2 <$ 


[-tan® , tan$ , -l] 
m m 


m 


m 


'■ 

1 t -9 


f- — 

A 

l 

V V 



-t 1 4> 


J 

V V 



0 -<t> 1 


6 

V V 


t . 


Substituting (2.23) and (2.24) in (2.21) and equating components, 


x = x - b 0 - 

v m m v 


■/l+tan 2 ® +tan 2 $ 


(tan® +0 -f tan4> ) 
m v v m 


m 


m 


y = b tj> + 
v mv 


m 


V 2 2 

l+tan © +tan 0 


(tan$ +4> -ty tan® ) 
m v v m 


m 


m 


h = b + 
v m 


m 


2 2 
l+tan © +tan $ 
m m 


(1-0 tan® — 0 tan<J> ) 
v m v m 


At nominal location and orientation 


x =y =0 , L =L , <J> =d =f =0 ; 

v v mo vvv 


thus 


$ =0 , x =L sin® , h =L co^3 +b 
o moo voom 

o 


( 2 . 22 ) 

for 

(2.23) 


(2.24) 

we get 

(2.25) 

(2.26) 
(2.27) 


(2.28) 
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Linearizing (2.25) - (2.27) about these nominal values gives 


x S - (b +L cos© )0 - L co9© 6 - sin© £ 

v mo ov o om om 

y £ (b +L cos© )4> + L cos© d> - L sin© f 

v mo ov o om o ov 

6h = cos© £ - L sin© 0 

v o m o o v 


where 

i =L -L , 0 =0 -0 and =$ -0 

mmo mmo mm 


(2.29) 

(2.30) 

(2.31) 


(2.32) 


Solving (2.29) - (2.31) for 0 , <l> and £ we have, 

mm m 


£ = sec © 6 h + L tan© 0 

m o v o o v 


0 = 

m 


sec © b 

y - (1+ — - sec © )ct> + tan© ty 

L v L ov ov 


o 
b 

0 = - sec 0 (sec © + —2)0 

m o o L v 

o 


sec 0 


— x + tan© 6 h 

L v o v 


(2.33) 

(2.34) 

(2.35) 


The angles between the cable carrying the hanging load and 
fuselage reference axis could be measured in the same way using 
elastically mounted ring. Thus 

V * V 9 v - -®v <1+ -P - T + T 

♦ = * -0 ■ -« (U ^ ^ 

£m £ v vi £ £ 


the 

another 


(2.36) 


The fuselage attitude could be obtained using an onboard vertical 

gyro. This gives measurements of 0 and 0 

v v 

All these measurements contain errors. We model them as additive 
white noise. The measurement vector z for the two systems can be 
written as 


z = Hx + v 

where the matrix H for the longitudinal and lateral systems is given 


15 



in Table II-2* , x is the state vector, and 

E£v(t)v T (T)] = R 6 (t-T ) (2.37) 


Table II-2 

Measurement Distribution Matrix H 


Longitudinal System 



Lateral System 



The measurements and the power spectral densities of the additive 
white noises are summarized in Table II-3. 


It is assumed that and h^ are determined accurately by 

separate measurements or another filter. 
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** This corresponds to a noise with RMS ,45m and correlation time .025 







2.6 System Definition Matrices for the Example Helicopter 

To illustrate the technique of designing the control system, we 
have taken a 6000 kg Sikorsky S-61 as the example helicopter. It 
carries a 2000 kg hanging load on a 20 m long cable, attached to a 
point 1,5 m below its center of mass. The measurement cable is 35 m 
long and is offset 25 m in the longitudinal direction from the refer- 
ence point on the ground. The state definition matrices F , G , 
and r are given in Table II-4 for the longitudinal-lateral motions 
of this system. Values of stability derivatives for this vehicle are 
taken from Hall and Bryson (H3). The drag coefficient on the hanging 
load is estimated assuming a reasonable size and drag coefficient and 
taking secant slope between zero and 7.0 m sec 1 wind velocity*. 

The control vector is redefined so that its two components are 

linear combinations of 0 and <t> as explained in section 2.4. If 

R R 

u = (0 +.03<t> , <f> -.030 ) the control distribution matrix is 

— R R R R 

modified from G to G as shown in Table II-4. 
o 

The decoupling approximation is made. The control u^ affects 
the longitudinal motion and the lateral motion. The state defi- 

nition matrices for the longitudinal and the lateral system are the 
appropriate portions of the matrices in Table II-4. Table II-5 gives 
H , the measurement distribution matrix. 

Figure 2.3 shows the eigenvalues of the uncontrolled longitudinal 
and lateral systems. Also shown are the corresponding eigenvalues from 
the coupled 12th order system. The eigenvalue identification is done 
using the eigenvectors. The eigenvalues from the coupled and the 
decoupled system are very near to each other showing that the decou- 
pling approximation is valid. 


The drag on a body depends on the square of the velocity. The slope 
of the drag vs. velocity curve at zero velocity is zero. 
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State Definition Matrices for One Version of Sikorsky S-61 Carrying a 2000 kg Hanging Load 
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Table I 1-5 


Measurement 

Definition 

Matrices 

for the 

Longitudinal System 

X 

T 

T 
H : 




e 

V 

-1.537 


-1.075 

1.0 

• 

e 

V 

0.0 


0.0 

0.0 

X 

V 

-.0348 


-.05 

0.0 

X 

V 

0.0 


0.0 

0.0 


0.0 


+.05 

0.0 

*J 

0.0 


0.0 

0.0 


z = [0 
— m 


\ 

0 ] 
V 

Measurement 

Definition 

Matrices 

for the 

Lateral System 

X 

i 

T 
H : 




4> 

V 

-1.052 


-1.075 

1.0 

D 

0 

V 

0.0 


0.0 

0.0 

y v 

.0348 


.05 

0.0 

y v 

0.0 


0.0 

0.0 


0.0 


-.05 

0.0 


0.0 


0.0 

0.0 


Z = [0 

— m 





20 




FIG. 2.3 Poles of the Open-Loop System 
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Chapter III 


DESIGN AND PERFORMANCE OF A PRECISION HOVER AUTOPILOT 
3.1 Introduction 

The control logic for the longitudinal-lateral motions is designed 
using quadratic synthesis and neglecting the cross-coupling of these two 
motions. The wind velocities are modeled as exponentially correlated, 
i.e., as first order continuous Markov processes. To test the decou- 
pling approximation, the control logic is used with the cross-coupled 
system model; the eigenvalues of this closed loop system are quite close 
to those predicted when neglecting the cross-coupling, though the eigen- 
vectors are strongly coupled. 

The performance of the system is first determined in the presence 
of winds assuming perfect state information (ideal case). 

Next, optimal filters are designed for the longitudinal and lateral 
systems using the measurements shown in Table I 1-3. The root-mean-square 
(RMS) values of state and control variables are found using the filters 
combined with the control logic designed before. The transfer functions 
of the filter-controllers (compensators) are determined and are shown to 
be non-minimum phase. 

If there is a steady wind disturbance the use of the control system 
above results in a steady state error. This steady state error can be 
avoided by using integral control, i.e., by adjoining a new state 
variable to the longitudinal system and another to the lateral system. 

For the longitudinal system the new state variable £ is defined as 

t = (3.1) 

This leads to zero error in the position of the hanging load in the 
presence of a steady wind. ! 
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3.2 


Controller Design for Exponentially Correlated Wind 


The correlation time of the wind velocity changes is usually 
comparable to or larger than the time constants associated with the 
controlled system. Therefore the wind disturbances cannot be modeled 
as "white noise". Hence, we model the head wind and crosswind as ex- 
ponentially correlated with a 5 sec time constant, i.e., 


V 

w 

X 



5 sec 


Therefore , 


V = - . 2V + n 
w w 'x 

X X 


(3.2) 


ti is white noise. If the RMS value of V is 7 m sec the 
‘x w 

X 

required spectral density of is 


„ 2X7 . 2 -3 

— = 19.6 m sec 

There is, of course, another equation like (3.2) for the crosswind 
disturbance. 


(3.3) 


If (3.2) is adjoined to the longitudinal dynamics equations the 
system is modified to 


_d_ 

dt 


X 


. 

— 

V 


w 


X 


— 



F ; r 


o ! -l/r 



~ X 


G - ' 


~CT 


~V 

w 

X 

4- 

0 

U + 

1 


_ — - 


_ ■— 


— -md 


(3.4) 


where x , u , F , r , and G are defined in Tables II-4 and II-5. 

We choose a control law such that u is a linear combination of x and 


V w , i.e., 
x 


u = Cx + c V 

- - W w 

X X 


(3.5) 


u is a scalar in this case; hence, C is a row vector and 
scalar. Substituting for u in (3.4) we have 


c 

w 

x 


is a 
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d 

dt 


(3.6) 



F+GC ! T+Gc 
"ix 


X 

+ 

0 

0 ! -l/T 

. c _ 


V 

w 

X _ 


i __ 


Notice that the closed loop eigenvalues depend on state definition 
matrices and the gains C 
choose gains C and c 


and do not depend on gain c . We can 

w x 


independently; C 


dynamics of the closed loop system, and c 
ance in the presence of wind. 


w„ 


to obtain the desired 
to achieve good perform- 


In the design procedure adopted here gains C are chosen using 

the quadratic synthesis procedure, i.e., by minimizing a performance 

index, J which is quadratic in state and control variables. The gain 

c can be chosen in several different ways. We have tried two of 
w x 

them here: 


(1) Using the augmented state model with the desired perform- 
ance index, quadratic synthesis gives the gain c in the 
same way as gains C . Quadratic synthesis strikes a com- 
promise between the amount of control required and state 
variable deviations. 

(2) In this procedure c is chosen to minimize the steady 

w x 

state RMS value of a desired linear combination of state 
variables (in this case the error in position of the hanging 
load). Thus, the control on wind can be tightened without 
making the closed-loop system too "fast”. This procedure is 
described in detail in Appendix C. 

The feedback gains on the state variables and wind velocity are shown in 
Table III-l for both longitudinal and lateral systems. The gains on wind 
velocity are computed using both procedures. 


The closed loop eigenvalues of the longitudinal and lateral systems 
are compared with the closed loop eigenvalues of the coupled longitudinal- 
lateral system (using gains from the decoupled systems) in Table 1 1 1- 1 . 

The eigenvalues are quite close showing that the decoupling approximation 
is reasonable. 


24 



Table III-l 


Controller for Exponentially Correlated Wind 



Feedback Gains 


c e 

V 

C 0 

V 

c 

! x 

V 

c. 

X 

V 

c 

x i 

c . 
X i 

C <D 

V 

C <t> 

V 

c 

y v 

c. 

y v 

c 

y i 

c . 

y f 

A. 

-1. 27 

277 

. 141 

. 157 

-.0296 

. 137 

-.834 

-.125 

-.137 

-.131 

.0249 

-.14 


Feedback Gains on Wind 



c 

c 

Technique 

w x 

w y 

1 

.00219 

-.00208 

2 

.00239 

-.00228 


Closed Loop Eigenvalues of Coupled and Decoupled Systems 


7th Order Longitudinal 
System 

[ sec '*'] 

7th Order Lateral 
System 

[sec 1 ] 

Complete 14th Order 
System 

[sec 

-.99 ± 1. 8 j 


-.93 ± 1 . 9 j 

-.58 ± l.Oj 


-.57 ± l.lj 

-.84 ± . 78 j 


-.91 ± . 62 j 


-2.2 ± 3 . 5 j 

-2.3 ± 3 . 5 j 


-.82 ± . 7 4 j 

-.72 ± . 7 9j 


-.53 ± l.Oj 

-.50 ± . 98 j 

-.2 

1 

-•2 

- . 2* 


* 

Double root 
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FIG. 3.1 Open-Loop and Close 
the Longitudinal Sy 
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FIG. 3.2 Open-Loop and Closed-Loop Poles of 
the Lateral System 
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Table III-2 


RMS Response with Perfect State Information 

(Wind RMS 7 m sec 1 longitudinal and lateral , 
correlation time 5 sec. ) 

Longitudinal System 


System 

Technique 


6 V 

x v 

x v 

x i 

x i 

u i 

Decoupled 

1 

. 00290 

.00495 

.0292 


asa 

.0169 

.00943 


2 

.00331 

.00601 

.0383 

|Q 9 

EH 

.0183 

.00948 

Coupled 

1 

.00329 

.00589 

.0381 

.0283 

.0367 

RES 

.0100 


2 

.00366 

.00682 

.0453 

.0324 

.0351 

m 

.0101 


Lateral System 


System 

Technique 


\ 

mm 

Uj| 


bb 

u 2 (^ R ) 

Decoupled 

1 

.00291 

.00697 

.0300 

.0251 

.0209 

.0141 

.00936 


2 

.00338 

.00873 

.0391 

. 0299 

.0183 

.0153 

.00939 

Coupled 

1 

.00335 

.00831 

.0389 

.0290 

.0346 

.0177 

.00989 


2 

. 00379 

.00984 



.0463 

.0334 

. 0330 

.0186 

.00992 


L Units in m , sec ; angles in rad] 
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Figure 3.1 shows how the poles of the longitudinal system move when 
the feedback gains of Table III-l are used. The poorly damped complex 
pole pair becomes better damped and faster. The unstable complex pole 
pair becomes stable and the two real poles change into a complex pair. 
The pole at -.2 corresponds to the wind model and cannot be moved by 
feedback since wind velocity is external to the system. The eigenvalue 
identification was done using the eigenvectors given in Appendix A. A 
similar graph for the lateral system is shown in Figure 3.2. 

The root-mean- square (RMS) response of the system is determined as- 
suming perfect state information (Bl, Chapter 15). Table III-2 shows 
the RMS values of the state and control variables when procedures (1) 
and (2) are used to find feedback gains on wind velocity. Also shown is 
the predicted performance when the gains from the decoupled systems are 
used with the coupled 14th order system. Procedure (2) gives a lower 
RMS value of the deviation in position of the hanging load though at the 
cost of higher RMS value of other state variables. The RMS errors are 
quite small for this strong gusty wind. 


3.3 Filter Design and RMS Response with Filter-Controller 

A filter is designed to estimate the state variables from noisy 
measurements of Table II-3. The sum of ground plus wind velocity can be 
measured by mounting velocity probes on the helicopter. However, it is 
not necessary to do so. The wind velocity can be estimated if all other 
state variables are observable from the measurements. We assume that 
velocity probes are not used and build a filter to estimate the state 
variables and the wind. If x is the estimated value of x , then the 
estimation equation for the longitudinal system is 


d 

dt 



f ' r 


" /s ^ 
x 


— 

G 

» 


1 

1 

lO 

+ 


0 . -l/T 


w 



_ i c^ 


L 


0 


K [z-Hx] 


(3.7) 


where z and H are defined for the longitudinal and lateral systems 
in equation (2.25) and Table II-5. The 7X3 gain matrix K is given 
in Table III-3. We define x , the error in estimate of x , thus 
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X 


X 


(3.8) 


- x 

The RMS estimation errors and the poles of the estimate error equations, 
are shown in Table III-3. 


Table III-3 

Filters for Longitudinal and Lateral System 


X 

longitudinal 


X 

^lateral 



’ 0 

V 

-.661 

-.482 

2.51 " 


0 

V 

-.670 

-.733 

3.37 ~ 

e 

V 

-1.61 

-1.07 

4.49 


0 

V 

-1.94 

-1.92 

7.65 

X 

V 

-17.8 

-.494 

-35.0 


y v 

18.6 

1.42 

25.0 

X 

V 

-.626 

2.60 

-26.7 


y v 

3.08 

-1.66 

21.6 

*! 

-16.9 

3.58 

-19.5 



16.6 

-3.68 

11.1 


-5. 13 

-. 142 

-9.25 


y i 

5.47 

.544 

6. 12 

V 

w 

X 

257 

203 

-717 

_ 


V 

W 

y 

-204 

-245 

796 

T 

z 

ii 

i — i 

CD 

3 

*711 

0 ] 
V 

T 

z 

= [<J> 

m 

V 

m 

<*> ] 
V 


Eigenvalues of the Estimate Error Equation (sec ) 


Longitudinal 

-1.1 + 

2.5j 

-.15 + . 82j 

-.34 + . 19j 

-2.3 

Lateral 

-1.7 + 

3.9j 

-.17 + . 81 j 

-.33 + . 23 j 

-3.3 


RMS Estimation Error 


e ! 

V 

\ , I 

0 ! 

V 

X 

V 

X 

V 

X i 

X 

£ 

V 

w 

X 

0 

V 

0 

V 

B 

a 

y i 

y i 

V 

w 

y 

.0062 

.0116 

.232 

.099 

.205 

.085 

3.59 

.0072 

.0246 

• 21 

.092 

. 187 

.083 

3.04 


[units in m , sec ; angles in radians] 
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FIG. 3.3 System Implementation 












Figure 3.3 shows how the system may be implemented. The complete 
system is divided into four parts: (a) controller, (b) system dynamics, 
(c) sensors, and (d) filter. 

The filter-controller can be regarded as a multi-input multi-output 
compensator. The transfer function matrices of the compensators and the 
pole-zero locations are given in Appendix B. The compensators for both 
longitudinal and lateral systems have a pole in the right half plane. 

Some elements of the transfer function matrices also have a zero in the 
right half plane. Hence, these are nonminimum phase compensators. 

The RMS response of the system is determined with the filter and the 
controller designed above in the presence of . 7 m sec 1 RMS longitudinal 
and lateral wind. Table III-4 summarizes the results for gain on wind 
from both techniques and also when the filter and controller gains from 
the decoupled systems are used on the coupled system. 

Table III-4 

RMS Response with Filter and Controller 
(Wind 7 m sec 1 RMS longitudinal and lateral) 

Longitudinal System 


System 

Technique 

6 

V 

6 

V 

X 

V 

X 

V 

X i 

*e 

v~~v 


1 

.0341 

.0736 

.403 

.313 

.470 

.241 

.0247 

Decoupled 

2 

.0342 

.0740 

.403 

.314 

.468 

.240 

.0249 


1 

.0355 

.0771 

.405 

.320 

.470 

.242 

.0253 

Coupled 

2 

.0357 

.0775 

.405 

.321 

.469 

.24 2 

.0255 

Lateral System 

System 

Technique 

V 

SR 1 ' a 

mm 

WM 

y i 


V*V 


1 

.0327 


.340 

. 288 

.394 

.202 

.0182 

Decoupled 

2 

.0328 


.339 

.287 

.393 

.202 

.0183 

Coupled 

1 

.0339 


.340 

.290 

.394 

.202 

.0187 


2 

.0341 


.340 

.291 

.393 

.202 

.0188 

i 

(Units in 

m , sec ; 

angles in radians) 
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There is considerable increase in RMS deviation of the hanging load 
position and control as compared to the case with perfect state informa- 
tion (see Table III-2). Nevertheless, the errors are still small for 
such a strong gusty wind. The coupled system gives nearly the same 
response. 

3 . 4 Integral Feedback 

The filter-controller of Sections 3.2 and 3.3 was designed assuming 
an exponentially correlated wind. In the presence of a steady wind, it 
produces a steady-state error. This steady-state error can be eliminated 
by the addition of integral feedback, which makes the position insensitive 
to steady disturbances or small changes in system parameters. 

We add a new state variable, £ » defined in equation (3.1) to the 
longitudinal system giving 



(3.9) 

(3.10) 

(3.11) 


The gains can be chosen either by quadratic synthesis or by a mix- 
ture of quadratic synthesis and pole assignment. If the closed loop 
system is stable and V w is a constant wind, then x can be brought 

X Xj 

to zero in the steady state, by proper choice of c^ . 

A filter is designed to estimate the state variable x . The new 
state variable, £ , can be estimated by using one of the two methods 
given here. 

(i) i = (3.12) 

A 

or (ii) £ = z = x. + v„ (3.13) 

x Z & 


With perfect state information the two methods are the same. 
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Table III-5 


Integral Control for the Longitudinal System 



(a) PSI* .00831 .0232 .0476 .0651 .0571 .0417 .233 .0122 

Filter 0.344 .0741 .401 .337 .481 .268 2.42 .0268 


(b) PSI .0179 0.0 .0398 0.0 0.0 0.0 1.78 .00926 

Filter .0378 .0704 .400 .330 .478 .264 2.99 .0256 

- - - 1 — - — ... i A 

2 -1 

(a) Wind with 19.6 m sec spectral density and small 
correlation time. 

(b) 7.0 m sec 1 RMS steady wind 
* Perfect State Information 
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Table I I 1-6 


Integral Control for the Lateral System 



C 4> 

V 

C; 

<t> 

V 

c 

y v 

c. 

y 

V 

c 

y i 

c. 

y i 

c 

n 

-.855 

-.125 

-.149 

-.138 

.0120 

-.169 

-.01 


Filter: 



4> 

V 

y v 

y v 

y i 


z 

-1.30 

-6.12 

18.6 

1.85 

16.8 

4.85 " 


n <b 

m 

-1.37 

-5.99 

1.43 

-2.82 

-3.46 

-.087 


4 >. 
im 

5.76 

23.7 

24.9 

26.2 

10.3 

8.53 


<t> 

V 


RMS State and Control 


Wind 

State 


0 

V 

4> 

V 

B 

y v 

y i 

y i 

1 


(a) 

PS I* 

.0127 

.0602 

.0492 

.078 

.0549 

.0411 

.0219 

.0128 


Filter 

.0352 

. 116 

.345 

.321 

.413 

. 233 

2.05 

.0220 

(b) 

PS I 

.0181 

0.0 

.0398 

0.0 

0.0 

0.0 

1.67 

.00924 


Filter 

.0373 

.0988 

.343 

.312 

.409 

.229 

2.64 

.0201 


2 -1 

(a) Wind with 19.6 m sec spectral density and small 
correlation time. 

(b) 70 m sec ^ RMS steady wind 


* Perfect State Information 
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Table II 1-5 shows the controller gains (obtained using quadratic 

synthesis) and filter gains for the longitudinal system. The filter 

2 -1 

was designed for a wind with 19.6 m sec spectral density and small 
correlation time. This table also shows the RMS state and control in 

2 _ i 

two different wind states: (a) wind with 19.6 m sec spectral density 
and small correlation time (approximated "white"), and (b) 7.0 m sec 1 
RMS steady wind. With filter the second method is used to estimate the 
augmented state variable. 

The response is excellent with perfect state information. In case 
(b) the root mean square (RMS) value of deviation in position of the 
hanging load is zero. With noisy measurements and a filter the RMS 
values of state and control variables increase but are quite reasonable. 

Similar data for the lateral systems are given in Table III-6. 

With integral control in the presence of steady wind disturbance, 
the wind velocity is not required for feedback and is not estimated. 

The above filter is designed to estimate other state variables. This 
filter is not optimal, since for a system in which the states are 
coupled through state transition, noise distribution or measurement 
distribution matrix the equations governing optimal state estimates 
are coupled. 

3 . 5 Summary 

An autopilot was designed for precision hover of a hanging load 
over a desired point. Satisfactory performance is obtained using a 
reasonable amount of control. 

The response is much better using the assumption of perfect state 
information (ideal case) than with noisy measurements and optimal filter. 
To improve performance significantly, it is necessary to have better mea- 
surements. Integral control will prove useful only if better measurements 
are available. 
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Chapter IV 


MOVING THE HANGING LOAD OVER SHORT DISTANCES 

4.1 Introduction 

In many applications of interest the load has to be carried over 
short distances. In such cases the time spent near the terminal 
points is large compared to travel time and it is advantageous to 
carry the load without retracting the cables even if it requires 
traveling only at moderate speeds. We assume that the helicopter 
travels sufficiently slowly that aerodynamic instabilities (coupled yaw 
and lateral pendulum motions) of the cable-load combination do not arise. 

In general the transfer is carried out in four stages: (a) zero 

velocity to cruise velocity at approximately constant acceleration, 

(b) cruise at constant velocity, (c) cruise velocity to zero velocity 
at constant deceleration, and (d) precision hover. For very short 
distances the cruise phase may be absent. 

In what follows we present a general scheme for tranferring the 
system from any state to hover at the origin of the coordinate system. 
Autopilot modes are chosen depending on the location of the system in 
state space. There are four possible modes: Acceleration, cruise, 

deceleration, and hover. At a switch curve the autopilot shifts from 
one mode to another and the system experiences a transient. We first 
discuss the mode selection logic and then the autopilot logic for each mode. 

4 . 2 Mode Selection 

The longitudinal system has six state variables. If one or more 

of the state variables 6 , & , 0„ , and 0„ has large initial 

v v Z i 

values, the system is placed in a stabilization mode. In this mode 
the position and velocity of the helicopter are not controlled but 
the pitch angle, the load cable inclination angle, and the corre- 
sponding angular rates are reduced to low values. It is necessary to 
do this for a simple switching diagram because acceleration or velocity 
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commands might otherwise produce unacceptable angular deviations for a 

short period of time. Thus as a starting condition we have small 

9 , 9 , 9 „ , and 9 . , and any x and x . 

v v i i v v 

The switching diagram is shown in Figure 4.1 in the x y , 

space. The cruise velocity is fixed at v and the commanded 

max 

acceleration at -a . Roughly, the reduced state space (x , x 

max v v 

space) is divided into four regions. The commanded mode is a function 

of x and x as follows: 
v v 

Acceleration Mode 
Deceleration Mode 
Cruise Mode 

Deceleration Mode 
Acceleration Mode 


If 


. 2 
X £ 
V 


2 

C £ 
v 


X > 

v 

X h 

t 

and 

and 

X 

< 

V 

V 

V 


max 

and 

X 

V 

> 

V 

max 

and 

X 

V 

- 

V 

max 

and 

X 

> 

-V 

V 

V 


max 

and 

X 

V 

> 

-V 

max 


If 




and 


1 x 5 x 

V 1 h 


Hover Mode. 


(4.1) 


The figure shows the details of the switching logic. It takes 
some time for the system to go into the new mode after mode switching 
is commanded. Therefore the system should be commanded to go into 
the new mode before it approaches the mode switching boundary of (4.1). 
The dotted lines in the diagram are the boundaries where the new mode 
is commanded. The switching curves for bang-bang type control have 
some dead space around them to avoid the problem of chattering. The 
dotted lines around the mode switching curves do not represent dead 
space alone. This is the space required for mode change to take 
place. 
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FIG. 4.1 Mode Selection Logic 





An example trajectory is shown. At the starting point, A , with 

the helicopter standing still, acceleration is commanded. The aircraft 

travels to point A' before constant acceleration is achieved. The 

Cruise Mode is commanded at B and the helicopter travels to point 

B' before the velocity is constant at v . This is followed by 

max 

the Deceleration Mode (C-C'-D) and the Hover Mode (D-D' ) . 


For the system under consideration the following parameters are 
chosen 


V 

max 

— 

15 m sec 

Av 

— 

2 m 

-1 

sec 

a 

max 

= 

1 m 

-2 

sec 

X h 

= 

1.0 

m 

V h 

= 

1.0 

m sec 


(4.2) 


4.3 Command Generating Model 


The command generating model produces commands to switch mode or 

change trim state with acceptable transient response. The desired 

changes are smoothed out so that the system can follow them with 

reasonable state variable deviations and control. The structure of 

the model depends upon the commands desired and the kind of transient 

response desired. In the model designed here, constant acceleration 

or velocity can be commanded. If a, is the desired acceleration, 

d 

the command velocity, v , and commanded acceleration, a , are given 

c c 

by 


v = a 
c c 


(4.3) 


a 

c 


1 

T 


(a -a ) 
d c 


(4.4) 
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Follower 


For velocity or acceleration commands, we are not interested in 

controlling position, so the longitudinal equations of 

motion can be written in terms of only five state variables, 0 , 

v 

9 , x , Q „ , and 0„ as shown in Appendix D. The controller can 

v v & £ 

follow a constant velocity or acceleration command in the "steady 
state ". In this "steady state", x (rate of change of acceleration), 
5 v and 9 are zero. Differentiating the state equations twice and 
then once gives 


and 

Now the 

a and 
c 


9 -9 = 0 

v £ 


(4.5) 


0~ 

__ _ 

.000254 

V 



0 . 


.000262 

L ll 

s 

_ 


( 4 . 6 ) 


state variables and the control variable are proportional to 


v^ in this steady-state : 


r el 

V 

= 

102 

-. 000254 " 

e 

V 


-.000254 

0.0 

X 

V 


0.0 

1.0 

CD 


-.102 

-.000262 

^ i 

•CD 

1 

S 

-.000262 

0.0 



= M I C 


( 4 . 7 ) 


u = -.0000158 a -.00132 v = Ny 
s c c — c 

The desired control law is 

u = u + C(x - x ) 
s — s 

The gains C are given in Appendix D. 


( 4 . 8 ) 


( 4 . 9 ) 
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FIG. 4.2 Implementation of Mode Switching Logic 












4.4 


Transient Response 


The transient response of the system was determined at mode switch- 
ing. The following four cases were considered with the time constant of 
the command generating filter equal to 2 sec: 

(a) Switching from hover to constant acceleration mode 
(Figure 4.3) 

(b) Switching from constant acceleration mode to constant 
velocity mode (Figure 4.4) 

(c) Switching from constant velocity mode to constant 
deceleration mode (Figure 4.5) 

(d) Switching from constant acceleration to constant 
deceleration mode (Figure 4.6) 

Good response was obtained in all cases as shown in figures (4.3) to 
(4.6). The closed loop frequency associated with pitch angle changes is 
higher than those associated with changes in inclination of the load 
cable or velocities of the vehicle and the hanging load. Thus after a 
mode is commanded, the pitch angle starts changing away from its previous 
trajectory first. The large transients in the rotor tilt in the begin- 
ning are due to fast changes in pitch angle and pitch rate. Notice that 

the initial transient response of the longitudinal tilt of the rotor NFP 

-2 

in figure 4.4 is similar to figure 4.5 (both involving a -1 m sec 

change in commanded acceleration) and is "negative" of figure 4.3 (in- 
-2 

volving a 1 m sec change in commanded acceleration). 

In all cases the maximum rotor tilt in the longitudinal direction 
(control) was less than 2° (the maximum control authority is about 12 s ), 
while the pitch angle never exceeded 8° (the nominal "steady state" pitch 
angle in acceleration or deceleration mode is about 5.7°). Switching 
from constant deceleration to constant acceleration produced the largest 
transient. 

4 . 5 Summary 

It is possible to move hanging loads over short distances safely and 
efficiently without retracting the cables. The cable-load configuration 
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FIG. 4.3 Transition from Hover to Constant 


Acceleration Mode 
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FIG. 4.3 (contd. ) Transition from Hover to Constant 


PITCH ANGLE, 

INCLINATION OF INCLINATION OF 

ROTOR NFP, 6 r (rad) LOAD CABLE (rad) VELOCITY (m sec ') 




FIG. 4.4 Transition from Constant Acceleration 
to Constant Velocity Mode 




PITCH ANGLE, 

INCLINATION OF INCLINATION OF 

ROTOR NFP,0 r (rad) LOAD CABLE (rad) VELOCITY (m sec 



TIME (sec) 


FIG. 4.5 Transition from Constant Velocity to 
Constant Deceleration Mode 
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INCLINATION OF 

ROTOR NFP, 0 R (rad) ANGLE (rad) VELOCITY (m sec 


HANGING 

LOAD 



FIG. 4.6 Transition from Constant Acceleration 
to Constant Deceleration Mode 






must be designed so that aerodynamic instabilities do not occur at 
velocities less than or equal to the cruise velocity. Studies by Gable 
and Wilson [ Gl] show that most practical hanging loads can be made aero- 
dynamically stable up to a reasonable speed. Some thoeretical studies 
[Pi] and experiments with one body shape and cable configuration [El] 
point to the contrary. However, these latter studies assume that the 
aircraft moves at constant speed and neglects the effect of the hang- 
ing load on the vehicle. This coupling is not negligible for modern 
crane helicopters and heavy lift helicopters, where the mass of the 
helicopter is comparable to that of the payload. 
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Chapter V 


PERFORMANCE UNDER OFF- DESIGN CONDITIONS 

5 . 1 Introduction 

In operation some of the parameters will differ from the values 
used in designing the control system, e.g., in off-loading container 
ships, all containers may not have the same mass (change in v) or loca- 
tion of the center of mass (change in £). In addition, the values of 
some of the system parameters may not be accurately known. Therefore 
it is prudent to study the performance of the system when the gains for 
the design condition are used with a different configuration. 

We have shown that the decoupling approximation is quite good for 
the purpose of' designing the controller and the filter and is reason- 
able for evaluating the performance. Also the behavior of the longi- 
tudinal and lateral systems is of the same character. Therefore it is 
sufficient to carry out a thorough investigation of the longitudinal 
system only. 

5 . 2 System Definition Matrices 

We believe the following four parameters are the most likely ones 
to change in normal operation: 

(a) mass of the vehicle, because of varying amounts of fuel 

m v +m^ 

(b) mass of the hanging load (or the ratio \ = — ) 

(c) the distance between the suspension point and the center 
of gravity of the hanging load 

(d) drag coefficient on the hanging load 

These parameters are changed one at a time and the matrices F , G , 

F , and H are computed. The mass of the vehicle is changed to 
5000 kg and 7000 kg (nominal value 6000 kg), the mass of the hanging 
load is changed to 1000 kg and 3000 kg (nominal value 2000 kg), the 
length of the cable, £ , is changed to 10 m and 30 m (nominal value 
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20 m) and the drag coefficient on the load is increased and decreased 
by a factor of ten from its nominal value of’ .0026 sec \ 

The open loop eigenvalues of the longitudinal system are shown in 
Table V-l. The character of the eigenvalues is the same. The mass of 
the hanging load and the length of the cable seem to be important 
parameters. 

5.3 Closed Loop System 

The nominal controller and filter gains obtained for the design 
configuration are used with the system matrices obtained for off- 
design conditions. Table V-2 shows the closed loop eigenvalues. The 
eigenvalues move but all the closed loop systems remain stable. The 
damping of all the roots is good except for one root pair for Y =1.1667 
(1000 kg hanging load). Filter eigenvalues are shown in Table V-3. 

All filters remain stable and the poles do not move much. It seems that 
it is possible to use the design filter and controller gains even with 
these wide fluctuations in parameters. 

Root mean square (RMS) state and control response in the presence 
of a 7m sec 1 RMS exponentially correlated wind (correlation time 
5 sec) are shown in Table V-4 with perfect state information (ideal 
case) and in Table V-5 with noisy measurements and filter with fixed 
gains. Except for the case with increased drag coefficient on the 
hanging load, the state and control responses do not deteriorate 
appreciably. With increased drag coefficient (last case) the RMS value 
of the deviation in the position of hanging load increases from .033 m 
to .481 m with perfect state information (PSI) and from .47 m to .77 m 
with noisy measurements and filter. We achieve suboptimal performance 
by using controller and filter gains which are optimal only under 
design conditions. The performance can be improved by choosing a 
different set of gains for each configuration but it does not seem 
necessary. In some instances this suboptimal performance is superior 
to the design case because those configurations are more favorable 
than the design configuration. 
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Open Loop Eigenvalues Under Off-Design Conditions 
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Eigenvalues of the Estimate Error Equation 
Off Design Parameters and Design Filter Gains] 
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Drag coefficient on hanging load 




Response with Perfect State Information 
Under Off-Design Conditions 
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** Design Condition 



Response with Noisy Measurements and Filter 
(Off-Design Conditions) 
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coefficient on hanging load 



5.4 


Summary 


Satisfactory system behavior is obtained by using a fixed set of 
controller and filter feedback gains even though some of the system 
parameters vary within wide ranges. It seems unnecessary to use 
different gains with nominal changes in system parameters. However, 
a new set of gains would be required with a completely new system 
configuration or measurement system. 
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Chapter VI 


SUMMARY 

A mathematical model is developed for the motions of a hovering 
or slowly moving helicopter carrying a hanging load. Autopilot logic 
is designed for position-hold (hover), velocity-hold, and acceleration- 
hold. 

It is shown that the longitudinal-lateral decoupling approximation 
is satisfactory for the purpose of autopilot design, except when a very 
tight control system is desired. 

The root mean square (RMS) response of the controlled system is 
quite satisfactory in the presence of steady and gusty longitudinal and 
lateral winds. Further improvements in hovering accuracy could be at- 
tained by using a better measurement system, since substantially better 
performance is obtained using the assumption of perfect state informa- 
tion (ideal case) under design conditions. 

The stability and performance of the system is examined for 
reasonable variations in system parameters. In no case does the system 
become unstable. The RMS response does not change appreciably except 
with a tenfold increase in drag coefficient on the hanging load. The 
behavior is satisfactory in all cases. 

It is possible to move the hanging load over short distances 
rapidly without retracting the cables. This involves controlling the 
oscillations of the load while moving at a low velocity. In case of 
failure of the automatic control system the pilot can still control 
the system. Some load shapes and cable arrangements are aerodynami- 
cally unstable (Gl, PI, El). These instabilities are caused by aero- 
dynamic forces on the hanging load and usually produce coupled oscil- 
lations involving yaw and lateral pendulum motions. Gabel and Wilson 
show experimentally that most hanging loads are aerodynamically stable 
over a wide speed range with a proper choice of cable configurations. 
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Before a load is moved, one should make sure that no aerodynamic 
instabilities exist up to the desired cruise speed. Some insta- 
bilities could probably be actively controlled if their nature were 
known in advance. However, it seems easier to avoid them in the 
speed range of interest. 

Further work is required into the nature of these aerodynamic 
instabilities. Then, autopilot logic to eliminate these instabilities 
could be developed, and, with adequate redundancy in the control 
system, it would then be possible to carry hanging loads at even 
higher speeds without retracting the cables. 
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Appendix A 


EIGENVALUES AND EIGENVECTORS OF THE OPEN-LOOP AND CLOSED-LOOP SYSTEMS 
A. 1 Open-Loop Systems 

Eigenvalues and eigenvectors of the open-loop coupled and decoupled 
systems are given in Table A-l in units of m , sec and radians . 


Table A-l 

OPEN-LOOP EIGENVALUES AND EIGENVECTORS 
(a) Longitudinal System 


Eigenvalue 

Eigenvector 

-1 

sec 

9 

V 

x 

V 


V 

-.11 ± l-7j 

.028 =F .llj 

-.040 T . 59j 

.025 + . 16 j 

0.0 


.0060 P . 0044 j 

.87 + .llj 

1.0 

0.0 

-.40 

-.11 

.75 

.55 ! 0.0 

0.0 

0.0 | -.71 

-.71 | 0.0 

-.20 

.0028 -.72 



(b) Lateral System 


Eigenvalue 

Eigenvector 

-1 

sec 

0 

V 

y v 


s 

-.67 + 2.8j 

-.095 + . 17 j 

-.078 =F . 33 j 

.013 + . 033 j 

0.0 


-.0065 + . 0044 j 


1.0 

0.0 

-.46 

.014 

.75 

.51 

0.0 

0.0 

0.0 

-.71 

-.71 

0.0 

-.20 

-.0028 

-.71 

-.66 

.14 


[Units m , sec ; angles in rad ] 
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Table A-l (Cont'd) 
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A. 2 


Closed-Loop Systems 


Table A-2 shows the eigenvalues and the corresponding eigenvectors 
of the decoupled and coupled longitudinal-lateral systems. 


Table A-2 

CLOSED- LOOP EIGENVALUES AND EIGENVECTORS 
(a) Longitudinal System 


Eigenvalue 

Eigenvector 

-1 

sec 

0 

V 

X 

V 


V w x 

-.99 + 1 . 8 j 

.11 T . 063j 

-.24 =F . 43 j 

.064 + . 020j 

0; 0 

-.84 + . 78j 

.066 =F . 028 j 

-.64 T . 60j 

.13 =F . 26 j 

0.0 

-.58 + l.Oj 

.058 =F . 04 Oj 

-.42 T . 74j 

.34 + . 076 j 

0.0 

-. 2 

.00026 

-.0043 

.0013 

1.0 


(b) Lateral System 


Eigenvalue 

Eigenvector 

-1 

sec 

0 

V 

y v 

y i 

H 

-2.2 + 3 . 5 j 

-.20 + . 097 j 

-.13 T „ 20 j 

.0074 + . 0087 j 

0.0 

- .82 + . 74 j 

-.063 + . 03 2j 

-.67 + . 61 j 

.12 =F . 29 j 

0.0 

-.53 + l.Oj 

-.050 + . 043 j 

-.40 =F . 77 j 

.39 T . 048j 

0.0 

-.2 

-.00026 

i 

-.0043 : .0013 

1 

1.0 


[ Units m , sec ; angles in rad ] 
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Appendix B 


TRANSFER FUNCTIONS OF COMPENSATORS FOR PRECISION 
HOVER IN THE PRESENCE OF CORRELATED WIND 

B. 1 Introduction 

As indicated in Chapter 3, a filter-controller can be considered to 
be a multi-input multi-output compensator. A possible implementation is 
shown in figure B. 1. In general this compensator has as many poles as 
the order of the system (except with pole-zero cancellations) and each 
term has at least one less zero than the number of poles. The compensa- 
tor transfer function matrix is 

A c = C{ S I - (F+GC-KH)} -1 K (B. 1) 

and is the transfer function between the measurements and control. 

B. 2 Compensators for Longitudinal and Lateral Systems 

The transfer function matrix of the compensator for both the 
longitudinal and the lateral systems is a 1X3 matrix. Each term has 
seven poles and six zeroes. The poles of the three terms are, of course, 
the same. The pole-zero location of the terms in the longitudinal system 
are shown in figures B. 2 to B.4 and of the lateral system are shown in 
figures B.5 to B.7. All elements in the transfer function matrices of 
the two systems have poles in the right half plane. Some elements have 
a pair of zeroes also in the right half plane. Such compensators are 
said to have non-minimum phase. 

It is clear that the closed-loop system has fourteen poles, seven of 
them can be identified with the controller and the remaining seven with 
the filter. All these fourteen poles lie in the left half plane and the 
closed-loop system is stable. 
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FIG. B.4 Compensator Transfer Function for 
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FIG. B.5 Compensator Transfer Function ^ for 
Lateral System 
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FIG. B.6 Compensator Transfe 
Lateral System 
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FIG. B.7 Compensator Transfer Function for 

Lateral System 
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These compensators cannot be tested directly since they are 
"unstable". One way to test them is to close the loop by using a 
simulation of the system. The response of the closed-loop system to 
initial errors and random disturbances can then be checked. 
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Appendix C 


A TECHNIQUE FOR FINDING FEEDBACK GAINS FOR 
CORRELATED DISTURBANCES 


C. 1 Introduction 

Consider a system 

x = Fx + Gu + Pw 

in which u is an m control vector and w is a p 
obeying the equation 

W = Ew + T) 

T) is white noise. If we choose 

u = Cx + C w 
— — w- 

in the steady state we have 

T T T 

fx + x f + r x +xr=o 

c xx xx c c xw xw c 

FX +XE T +PX =0 

c xw xw c ww 


(C.l) 


disturbance vector 


(C. 2) 


(C. 3 ) 


(C. 4 ) 
(C. 5 ) 


where 


3fe 

With noisy measurements and optimal filter (See (C.6) for definition) 

X = X^S/N + X 

xx xx xx 

where x is the estimate of x and x is the estimation error. X 

— • — — xx 

depends on state and measurement noise. Choice of gains C^ only 

changes X^~ . The equation governing x is 
xx 

x = Fx + Gu + H^ - KH~ 

x 
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(C. 6 ) 


X= E(ab T ) 
ab 


F = F + GC 
c 

r = r + gc 

c w 


We choose gains C first to obtain the desired closed-loop 

dynamics. The gains can subsequently be chosen to reduce the 

effect of wind as much as possible, i.e., if Tx are q independent 

linear combinations of the states we are interested in, C could be 

w 

chosen to minimize 


rp 

L = Tr E{ (Tx) DTx} s.s. 


= Tr{ T DTX } 
xx 


(C. 7 ) 


where D is qXq positive definite weighting matrix. Thus the 
Hamiltonian is 


H = L + Tr 


{ f x + x f t + r x T + x r T }A 

C XX XX c c xw xw c 


+ {FX + X E T + rX }x 

C XW XW r WW 


(C. 8 ) 


Here A and X are Lagrange multiplier matrices. The required opti- 
mization equations are 


(a) (C. 4 ) and (C.5) 

^ U T T 

(p) = 0 => F c A + AF + T DT = 0 

xx c 

d H T 

5 — = 0 2 AT + F X + A.E = 0 

o X c c 

xw 


^ u T T T 

(V ) * 7 ; =0=^2GAX + G A X =0 

o C xw ww 

w 


(C. 9) 


(C.10) 


(C.ll) 


Comments on solving the optimization equations: 

1. Equation (C.9) can be solved for A . 

2. Equation (C.4) can be solved in the end since no other equation 
involves X 

xx 
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3. (C.5), (C. 10) and (C. 11) are (2n+m)p equations in (2n+p)m 

unknowns X , \ and C 
xw w 

4. These equations assure a solution for well posed problans (i.e., 

negative definite F and negative semidef inite E , positive 

c 

definite X and finite, nonzero G ). However, the solution 

ww * ’ 

may not be unique. Any nonuniqueness can be resolved by minimizing 

T„ ’ 

u Bu . 

5. If the space spanned by vectors in G contains the space spanned 

by vectors in P , then one solution is to choose C such that 

w 

r + GC is zero, in which case L = 0 . 
w ’ 


C. 2 Random Bias Disturbance 

For bias disturbances the matrix E is zero and T] = 0 . Then 

since F is negative definite and hence invertible, equations (C.5) 
c 

and (C. 10) can be solved directly. 


x = -f 1 r x 

XW C C WW 


-1 T 

X = -2F A AT 
c c 


(C. 12 ) 


Substituting (C.12) and (C.9) in (C.ll) gives 


MC + N = 0 
w 

where 

T -T T -1 
M = G F T DTF G 
c c 

T -T T -1 
N = G F T DTF T 
c c 


(C. 13 ) 


(C. 14 ) 


The solution is unique if M has rank m , i.e., if 

T -T T 

rank(M) = rank(G F T ) = m (C.15) 

c 


The rank of M cannot exceed m , which is the rank of G . If there 
are exactly m disturbances and M has rank m 


-1 -1 

C = -(TF G) TF T 
w c c 


(C. 16) 
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C. 3 


Independent Disturbances 


The design procedure is simplified when the p disturbances are 

independent, i.e. , E and are diagonal (this generalizes to the 

case where E and X are block diagonal with the same block sizes). 

ww 

Then the gains on each disturbance can be chosen independently. To 
find gains on any particular disturbance we proceed as follows. 

In our general formulas, we consider p equal to one and E and 

X scalars. (C.5) and (C.10) become 
ww 

(F +EI )X + T X = 0 
c xw c ww 

(F +EI) T \ + 2AT = 0 
c c 


or 


X = -(F +EI) 1 r X 
xw c c ww 


X = -2(F +EI)“ T AT 
c c 


(C.17) 


Substituting in (C. 11) we have 

MC + N = 0 
w 

where 

M = G T F T (T T DT - 2EA)F G 
ce ce 

N = G T F T (T T DT - 2EA)F T 
ce ce 

and F = (F + El)" 1 
ce c 

Nonuniqueness can be resolved as in section (C.l). 


(C. 18) 


(C. 19) 
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Appendix D 


VELOCITY HOLD AUTOPILOTS FOR LOW SPEED 


D. 1 Design 

A velocity-hold autopilot maintains the vehicle speed near a desired 
value. With no constraint on time, we are not interested in errors in 
horizontal position. Thus the state equations of Chapter II can be re- 

. # . x 

written in terms of the state variable vectors, (0 ,6 ,x ,0„,0.) 

T v’ v* v’ £’ £, 

and (<t> ,4> ,y .$„,<!>) . Now in addition to the measurements of four 

V V v i i 

angles and <t> v we need measurements of longitudinal and 

lateral velocities. Noise in measurement of velocity is assumed to be 

2 -1 

white with power spectral density .005 m sec . The state definition 
matrices and the measurement distribution matrices are shown in Table 
D-l. We estimate that in the speed range of interest the drag coefficient 
terms on the vehicle and the load do not change substantially. Small 
changes will not deteriorate the performance very much as has been shown 
for the precision hover autopilot. 

The lateral velocity is usually maintained at zero while the 
longitudinal velocity is held at the desired cruise speed. 

Both longitudinal and lateral systems are augmented with a state 
variable corresponding to wind velocity. The control systems are 
designed using quadratic synthesis. The feedback controller gains 
are presented in Table D-2. Also shown are the open-loop and closed- 
loop eigenvalues. All modes of the closed-loop system have adequate 
damping. 

Filters are designed based on the measurements given above in the 
presence of wind with RMS intensity 7.0m sec * and correlation time 
5 sec . Filter gains and characteristic roots of the estimate error 
equation are given in Table D-2. 
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LONGITUDINAL SYSTEM 
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(Units in m , sec ; angles in rad 






Table D-2 


Controllers and Filters for Velocity Hold Autopilots 


long. 

= { 

100 

e 2 

V 

.2 

+ x + 

V 

200 


+ 100 

2 

U 1 

^ s. s. 

lat. 

= £ 

100 

♦: 

.2 

+ y + 

V 

200 

^ 2 

+ 100 

2 

U 2 

^ s. s. 


Feedback Gains 


o 

CD 

< 

> 

•CD 

U 

c. 

X 

V 


c e 

C 

w 

X 

C $> 

V 

El 

C. 

y v 

S 


C 

w 

y 

-1.70 

-.431 

.0981 

-.714 

-.901 

.00174 

-1.39 

-.242 

-.0981 

-.912 

-.824 

-.00172 



.0192 

-.0166 

-.872 

.260 

.0392 

29.0 

K? 

3.19 

5.68 

-19.6 

.0770 

.541 

-886 

long. 

-.0600 

-.0914 

.908 

-.0107 

-.0315 

20.4 


-.0190 

-.0802 

.793 

.263 

.0391 

-45.2 

k" = 

4.11 

8.85 

16.1 

-.076 

.145 

985 

lat. 







.0493 

.0580 

.834 

.0097 

.0302 

11.8 


Eigenvalues (sec ) 


Longitudinal 

System 

Open- Loop 

-.11 ± 1. 7 j 

m&mm 

-.41 

-.20 

Closed-Loop 


-.40 + . 64 j 


-.20 

Estimate Error 

-.95 + 2. 3 j 

-.25 + . 79j 

-1.9 

-.70 

Lateral 

System 

Open- Loop 

-.67 + 2. 8j 

.098 + . 27 j 

-.46 

-.20 

Closed-Loop 

-.34 + 4. 4 j 

-.38 + . 63 j 

-1.45 

-.20 

Estimate Error 

-1.5 + 3 .7 j 

-.24 + .77 j 

-2.9 

-.67 


(Units in m , sec ; angles in rad ) 
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D. 2 


Performance 


The performance of the system is computed with perfect state 
information (ideal case) and with noisy measurements and filter. Table 
D-3 shows the results. The performance is satisfactory. Major improve- 
ments in RMS response can be affected only by improving the measurements. 


Table D-3 

RMS Errors for the Longitudinal 
and Lateral Systems 



e 

V 

b 

V 

X 

V 

a , 


u 

1 

Perfect State 
Information 

.00173 

.00148 

.0171 

.00197 

. 00085 

.00931 

Noisy Measurement 
and Filter 

.0143 

.037 9 

. 136 

.0095 

.0094 

.0189 



<t> 

V 

mm 

^v 


0 

*J 

U 2 

Perfect State 
Information 

.00173 

.00217 

.0166 

.00194 

. 00085 

.00928 

Noisy Measurement 
and Filter 

.0140 

.0566 

.119 

.0081 

.0093 

.0146 


(Units in m , sec ; angles in rad ) 


81 
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